Suppression of the Shear Viscosity in a "semi" Quark Gluon Plasma 
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We consider QCD at temperatures T near T c , where the theory deconfines. We distinguish 
between a "complete" Quark Gluon Plasma (QGP), where the vacuum expectation value of the 
renormalized Polyakov loop is near unity, essentially constant with T, and the "semi"-QGP, where 
the loop changes strongly with T. Lattice simulations indicate that in QCD, there is a semi-QGP 
from below T c to a few times that. Using a semi-classical model, we compute the shear viscosity, 
r], to leading order in perturbation theory. We find that near T c , where the expectation value of 
the Polyakov loop is small, that rj/T 3 is suppressed by two powers of the loop. For heavy ions, this 
suggests that during the initial stages of the collision, hydrodynamic behavior at the LHC will be 
characterized by a shear viscosity which is significantly larger than that at RH1C. 



PACS numbers: 

The collisions of heavy ions at RHIC have demon- 
strated clear signals for a novel regime. Much interest 
has focused on collective properties, especially elliptical 
flow, which appear to be well described by a system in 
which the (dimensionless) ratio of the shear viscosity, 77, 
to the entropy, s, is small [1, 2, 3, 4, 5, 6, 7, 8, 9, 10]. 
These results do not agree with the expectations of a 
weakly coupled plasma, and have been described as a 
"strong" Quark Gluon Plasma (QGP) [3, 10]. 

It is reasonable to expect that the running coupling in 
QCD, a s (T) = g 2 (T) / (An) , might be large near the criti- 
cal temperature, T c . Typically, one expects the coupling 
to be in a non-perturbative regime for momenta less than 
~ 1 GeV. Since the transition temperature is ~ 200 MeV 
[11, 12, 13, 14], the coupling could well be large at several 
(and maybe many) times T c . 

Indeed, perhaps the coupling is so large that the rel- 
evant limit is of infinite coupling. It is possible to com- 
pute in a J\f = 4 supersymmetric SU(iV c ) theory when 
the number of colors, N c , and a s N c , are both infinite 
[15, 16, 17]. The Af = 4 supersymmetric theory is con- 
formally invariant, so rj/s is independent of temperature, 
and is small, = 1/4-7T [15, 16]. If QCD is analogous to 
the Af = 4 theory [17], then for some region above T c , 
77/ s should remain small and not change markedly with 
temperature. 

To have a coupling which is relatively moderate in 
strength at T c would be exceptional. At least for the 
pressure, this might occur because of the ubiquitous fac- 
tors of 2ir which accompany the temperature T in the 
imaginary time formalism [18]. Using an effective theory 
in three dimensions [19], a two loop computation shows 
that this does, in fact, occur: at T c , the effective coupling 
is only af w 0.3 [20]. 

The challenge is then to understand how the confining 
transition, with a large decrease in pressure, occurs for 



moderate coupling [21]. We find it useful to view decon- 
finement as the ionization of color charge. Without dy- 
namical quarks, in the confined phase there is no ioniza- 
tion of color. Conversely, far into the deconfined phase, 
color is completely ionized, either with or without quarks. 
In a non-Abelian gauge theory, the expectation value of 
the (renormalized) Polyakov loop characterizes the de- 
gree to which color charge is ionized. Lattice simulations 
[11] find that the Polyakov loop is small near T c , and near 
one at a few times T c [12, 13]. This regime, which we view 
as one of partial ionization, coincides with the drop in the 
pressure (relative to the ideal gas term). For want of a 
better term, we refer to this as the "semi" -QGP. Above 
a few times T c , there is a "complete" QGP, where both 
the renormalized Polyakov loop, and the pressure/T 4 , are 
essentially constant. 

In this paper we consider how the shear viscosity 
changes in the semi-QGP. Wc make numerous drastic 
assumptions. If the coupling is moderate even down to 
T c , perhaps we can treat the semi-QGP by means of a 
semi-classical approximation. Our ansatz is extremely 
simple, just a background field with constant Ao, Eq. (1). 
We then compute to leading order in a s and log(l/a s ) 
[22, 23], for an infinite number of colors and flavors [24]. 
We ignore the change of 77 with a s , to concentrate on 
how it changes as the expectation value of the Polyakov 
loop decreases. We find that when the loop is small, that 
the shear viscosity is suppressed by two powers of the 
Polyakov loop. This implies that 77/T 3 decreases signifi- 
cantly in the semi-QGP as the theory cools, from a few 
times T c down to T c . Since rj ~ 1/c^ at small a s , Eq. 
(3), including the running of the QCD coupling could re- 
inforce this trend. Such a large decrease near T c is very 
different from the strong QGP, where 77/s changes little 
[3, 10], if at all [17], with temperature. 

There are numerous examples of non-relativistic sys- 
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terns which exhibit a minimum in the shear viscosity near 
T c [8, 15]. The present analysis suggests how this might 
arise dynamically in QCD; see, e.g. [3]. Our result is also 
reminiscent of the "anomalous" viscosity of an Abelian 
plasma, where the viscosity changes because of a back- 
ground field, although our specific background field is not 
like that of Abelian plasmas [5] . 

The Polyakov loop represents the propagator of an in- 
finitely massive test quark, which we assume is in the 
fundamental representation [12, 13, 21, 25, 26, 27, 28]. 
Its magnitude can be viewed as the probability for the 
test quark to propagate. This is near one at asymptoti- 
cally high temperature, where the plasma is completely 
ionized. 

In a SU(iV c ) gauge theory without dynamical quarks, 
in the confined phase the propagator for a test quark 
vanishes identically. This is because the Polyakov loop 
carries Z(N C ) charge, and the confined phase is Z(N C ) 
symmetric. Thus without dynamical quarks, there is 
absolutely no ionization of Z(N C ) charge below T c , only 
above. 

The expectation value of the renormalized Polyakov 
loop, (£), is extracted from that of the bare loop after 
a type of mass renormalization [12]. For a SU(3) gauge 
theory without quarks, from Fig. (1) of Rcf. [12], (I) = 
below T c ; at T+, (£) ~ 0.5; it crosses one at ~ 2.8 T c , and 
reaches ~ 1.1 by ~ 4T C . This value is then constant from 
~ 4T C up to to highest temperature measured, ~ 12 T c . 
(In perturbation theory, (£) exceeds unity.) Thus in a 
SU(3) gauge theory without quarks, there is a confined 
phase below T c , a semi-QGP from (exactly) T c + to about 
~ 4T C , and a complete QGP above that. By its nature, 
the boundary between the semi- and the complete QGP 
is not precise. 

In QCD the Polyakov loop is no longer a strict or- 
der parameter, since dynamical quarks also carry Z(N C ) 
charge. Thus even below T c , the Z(N C ) charge of a test 
quark is shielded by the thermal ionization of dynamical 
quark anti-quark pairs. With sufficiently many flavors of 
quarks, this could happen even at rather low tempera- 
ture. 

At present, however, numerical simulations indicate 
only a modest ionization of color below T c , at least for 
three colors and 2 + 1 flavors of dynamical quarks. From 
Fig. llofRef. [13], the expectation value of the Polyakov 
loop is very small below 0.8 T c ; it then rises to ~ 0.3 at 
T c , and is near one at ~ 2 T c . We take this to show that 
in QCD that there is a semi-QGP between ~ .8T C to 
perhaps ~ 2 — 3T C . Thus with dynamical quarks, there 
is a semi-QGP in both the hadronic phase, from ~ 0.8 T c 
to T c , and in the deconfined phase, from T c to ~ 2 — 3 T c ; 
again, there is a complete QGP above that. 

We characterize the semi-QGP by the following ap- 
proximation [21, 25, 26, 27, 28]. The Polyakov loop is 
the trace of a straight Wilson line in imaginary time. An 
expectation value for the Polyakov loop which is not near 



one implies a non-trivial distribution for the eigenvalues 
of this Wilson line. We thus expand about a constant, 
background field for the timelike component of the vector 
potential, 

At = Q/g . (1) 

where g is the coupling constant for an SXJ(N C ) gauge 
theory. The matrix Q is diagonal in color space, with 
the indices those for the fundamental representation, a = 
l...N c . The Wilson line is L = exp(zQ/T), and the bare 
Polyakov loop is I = tr L/JV C . Although we expand about 
a given Q, what is physically relevant is the distribution 
of Q's. The moments of this distribution are given by 
powers of the Wilson line, trL 3 /N c , j = 1 .. . (iV c — 1). 

In a complete theory we would sclf-consistently deter- 
mine the effective Lagrangian which gives the distribu- 
tion of the Q's- Notably, this would enable one to com- 
pute the pressure (and thus the entropy) as a function 
of the Q's; e.g., versus (£), etc. Some results for the Q- 
distribution can be obtained analytically when the spatial 
volume is a small sphere [26]. In a large volume, its form 
is unknown [21, 27]; its determination through numerical 
simulations of effective theories appears promising [28]. 
In lieu of this, we take the first moment of the Wilson line 
from lattice simulations in the full theory [12, 13], and 
choose two different forms for the complete eigenvalue 
distribution. We find, unexpectedly, that at least the 
shear viscosity is rather insensitive to which distribution 
we take. 

The shear viscosity is computed using a Boltzmann 
equation [22, 23] in the presence of the background field, 
Q. In the imaginary time formalism, the Euclidean four 
momenta are = (po,p), where po is an even (odd) 
multiple of ttT for bosons (fcrmions). To compute with 
Q 7^ we use the double line notation of 't Hooft [24]. 
To simplify the calculation we take N c and Nf to be infi- 
nite, although the method can be used at finite values 
[29]. How the Q's enter depends upon the represen- 
tation of the gauge group. Quarks in the fundamental 
representation have one color line, so the four momenta 
has one color index, = (po + Q a ,p), where Q a is 
the diagonal component of the matrix Q (we also use 
P^ a = p n — Q a ). Adjoint gluons carry two color lines, and 
their momenta has two color indices, P^ b = {po + Q ab , p), 
defining Q ab = Q a — Q b . Each component Q a is typically 
a non-integral multiple of 2irT, and so Q a ^ cannot be 
shifted away without violating the appropriate boundary 
conditions in imaginary time [22]. Q a is like an imagi- 
nary chemical potential for color charge, in the space of 
diagonal (mutually commuting) generators. 

Amplitudes are analytically continued from imaginary 
to real time in the usual fashion, except that now it is 
necessary to analytically continue using energies which 
carry color, po + Q a — > — ioj a for quarks, and po + Q ab — > 
— iu} ab for gluons [30]. The mass shells remain on the 
light cone, to = ±E, E = \fp*. 



3 



For quarks, what typically enters are distribution func- 
tions n(E — iQ a ), where n(E) is the usual Fermi-Dirac 
statistical distribution function. To compute, expand 
this as [26] 

e (E-iQ«)/T + 1 "2^ > ■ W 

3=1 

The first term, ~ exp(— E/T), represents the Boltzmann 
approximation to the quantum distribution function, and 
is accompanied by exp(iQ a /T). Expansion to j th order 
brings in a factor of ~ exp(— jE/T), and is accompanied 
by exp(ijQ a /T). For a given process, the moments of 
the Wilson line which enter depend upon the detailed 
routing of the color indices and the like. As an example, 
consider the trace of the quark propagator: then the first, 
Boltzmann term involves the trace of exp(iQ a /T), which 
is the Polyakov loop, £; terms to j th order involve the j th 
moment of the Wilson line, trlA 

The gluon propagator is similar, except that n(E — 
iQ ab ), enters, where n(E) is the Bose-Einstein statisti- 
cal distribution function. Again, as an example consider 
summing over the indices of the gluon propagator. To 
avoid taking the trace, which is part of the correction in 
1/N C , one sums separately over a and b. The first, Boltz- 
mann term involves the traces of cxp(iQ ab /T), which be- 
comes | tr L | 2 ; terms to j th order become |tr L- 7 1 2 . 

In computing pcrturbatively about a trivial vacuum, 
one naturally divides the momenta into hard momenta, 
whose components are ~ T, or soft momenta, where both 
the energy, u>, and spatial momenta, p, are ~ gT [22]. 
This remains valid at Q a ^ for the colored energies, oj a 
and Lu ab ; the Q's themselves are ~ T, and so hard. 

We have computed the hard thermal loops (HTLs) in 
the quark and gluon self energies for Q ^ [29]. HTL's 
are the dominant contributions when the external mo- 
menta are soft. To leading order, ~ g 2 7 the dominant 
term involves an integral over hard momenta, times an 
angular integral over soft momenta [22] . The dependence 
on the Q's only affect the integral over the hard mo- 
menta, through the change in the statistical distribution 
functions; the angular integral is unchanged. 

For the quark self energy, the only change in the HTL is 
the change in the Dcbye mass, which now depends upon 
the direction in color space. For the gluon self energy, 
there is the usual HTL, again modified by the change 
in the Debye mass. Besides the HTL, which is ~ g 2 T 2 , 
there are also terms ~ g 2 T 3 Q/uj. These new terms are 
directly proportional to tadpole terms, and arise because 
Q ^ induces a net color charge. To understand the 
terms ~ T 3 , consider a Z(N C ) interface, which is modeled 
by a spatially dependent Q [25]. There such terms are 
physical, because the interface induces a nonzero color 
electric field [29]. For the present problem, however, we 
assume that there is no net color charge in the vacuum, so 
the only change in the soft gluon propagator is the change 




FIG. 1: The function /(L) of Eq. (3), versus L 'Step' and 
'GW denote Q distributions with a simple step-function and 
that in the Gross- Witten matrix model [26, 27], respectively. 



of the Debye mass. (If the gluon HTL does include a 
term ~ T 3 , then for processes involving exchange of a 
soft gluon, in Eq. (3) log(l/g) -> log(l/ 5 2 / 3 )). 

We now outline the calculation of the shear viscosity 
[29]. To leading order, both in g 2 and in log(l/g), we 
find that the result can be written as [22, 23, 29] 

W = 4 i C n / , /(L) • ( 3 ) 
T g 4 log(l/.g) 

The constant c v depends upon the number of colors, N c , 
and flavors, Nf [23]. As noted before, for calculational 
reasons we compute for N c = oo, with Nf/N c fixed, im- 
plicitly assuming that /(L) is relatively insensitive to N c . 

At this order, the viscosity is determined by the scat- 
tering of 2 — > 2 particles, where all particles have hard 
momenta. They interact through the exchange of a sin- 
gle, soft field in the i-channcl. For the soft field, we use 
the HTL approximation. The Debye mass changes when 
L / 1, but this doesn't enter at leading logarithmic order. 

For the pure glue theory at infinite N c , to leading or- 
der in g 2 and log(l/g>) [22, 23], two scattering processes 
contribute to the collision term: P? b + P| c -> P 2 ad + Pf c 
and Pf + Pf -> P£ d + Pf. This is nearly forward scat- 
tering, with the spatial momenta pi ~ P2, and p% s=s p^. 
The momenta of the exchanged gluon, p\ — P2 = Pa ~P3, 
must be soft to give log(l/g). 

Under these approximations, the dependence upon the 
background field Q only enters through the integrals over 
the statistical distribution functions for the hard fields, 
Pi and P2- While the statistical distribution functions 
are complex when Q ^ 0, after summing over both emis- 
sion and absorption processes, all contributions to the 
stress energy tensor are real. For the shear viscosity, the 
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integrals which enter are 

P OO 

/ dpp i n(p-iQ ab ) , 
Jo 

/>oo 

/ dp p> n(p - iQ ab ) (1 + n{p - iQ ad )) , (4) 
Jo 

where j = 2 and 4. One then expands the distribution 
functions as in Eq. (2), to obtain a power series in mo- 
ments of the Wilson line. With quarks, there are more 
diagrams, as the distribution functions between quarks 
and gluons mix [23]. The integrals are like Eq. (4), but 
also involve n(E — iQ a ). (We also take xip) ~ P 2 [23]: 
this is valid to to < 1.0% for Q = 0, and is exact at small 
£, where Boltzmann statistics applies.) 

As written in Eq. (3), the result is a function of the 
Wilson line, /(L), times the usual pcrturbative result; 
thus / = 1 when L = 1. This simple form is not valid 
beyond leading logarithmic order. If I = tr L/iV c is small, 
and dominates higher moments, we find the result van- 
ishes like the square of the loop: 

/(L)« 02(A) ^±^W *<1, (5) 

where A = N f /N c : o 2 (0) « 3.31, a 2 (l) « 1.01. 

The behavior at small £ can be understood as follows. 
In the pure glue theory, when gluons with momenta Pf 6 
and P$ c scatter, the Q b charges cancel, so summation 
over a and c gives a collision term ~ g A £ 2 (times log(l/g), 
which we suppress for brevity). There is no such cancel- 
lation for the scattering of gluons with momenta Pf 6 and 
Pg d , which is only ~ g A £ A ■ The source term [23] is like 
the trace of the gluon propagator, ~ £ 2 , and so without 
quarks, at small I the gluon contribution to the shear 
viscosity is %luon ~ (£ 2 ) 2 /(g A £ 2 ) ~ £ 2 /g A . 

With quarks, the collision term is dominated by a 
quark scattering off of an anti-quark, Pf + P 3 ~ a — > P 2 + 
P 4 _c . The Q a charges cancel, so this is ~ g A £° at small £. 
If a quark scatters off of a quark, Pf + P| -> P| + Pf, 
the Q charges don't cancel, and the scattering is ~ g A £ 2 ■ 
Gluons also scatter off of quarks, by exchanging a gluon: 
Pf 6 + P| -> P 2 ac + Pf; the Q b charges cancel, so this is ~ 
g A £. The mixing of the gluon and quark distribution func- 
tions [23] can be neglected at small £: quark anti-quark 

annhilation, Pf + P 3 ~ 6 -> P 2 ac + P 4 c( ~ 6) , is - g A £ 2 , while 
Compton scattering, Pf +P| C -> P 2 ab + Pf, is ~ £ 4 £ 3 . For 
the source term, the quark contribution is like summing 
over the quark propagator, ~ £, so the quarks contribute 
to the viscosity as 77q U ark ~ (£) 2 / '(g 4 1 £°) ~ £ 2 jg A . In the 
presence of quarks, gluons contribute to the viscosity as 
%iuon ~ (£ 2 ) 2 /(g A £) ~ £ 3 /g A , down by - £ to the quark 
contribution. 

To obtain results valid for all values of £, some assump- 
tion about higher moments must be made. We used two 
forms. The first is to take a simple step function, of width 
P, about the origin, so that trL n /iV c = s'm(n(3)/(n/3). 



The other is to take a Q distribution as in the Gross- 
Witten matrix model [26, 27]. By the methods described 
above, it is straightforward to obtain results, although 
their analytic form is unwieldy These forms can be eas- 
ily evaluated numerically, though, as shown in Fig.(l). 
We find that /(L) is insensitive to the assumption about 
higher moments, changing by at most a few percent over 
the entire range of £. We find that with both eigenvalue 
distributions, that the collision term generates a cusp in 
/(L) near I = 1. We expect that this non-analytic be- 
havior will be washed out by corrections to higher order, 
which enter for Q ~ gT. 

We conclude with some general comments. The 
Polyakov loop is proportional to the propagator of an 
infinitely massive test quark, and as such, has no direct 
relation to the propagation of dynamical fields. In our 
ansatz, however, the quasiparticles are fluctuations about 
the background field in Eq. (1). When the expecta- 
tion value of the Polyakov loop is small, this background 
field universally suppresses the propagation of any col- 
ored field. For heavy fields, this is reasonable, but for 
light fields, it is nontrivial. If the light fields have hard 
momenta p >T, then suppression by the Polyakov loop 
is the dominant effect, with other corrections down by 
powers of g. At soft momenta, p ~ gT, light fields are 
not only suppressed by a (small) Polyakov loop, but al- 
tered by the change in their hard thermal loops. 

In general, the shear viscosity is proportional to the ra- 
tio of a source and a collision term [22, 23] . In the absence 
of a background field, the source term is of order one, and 
so the shear viscosity can only be small if the coupling 
constant, and so the collision term, are large. This is the 
central idea which motivates the strong QGP [3, 10]. In 
contrast, in our analysis of the semi-QGP, the quasipar- 
ticles are weakly coupled, but are fluctuations about a 
nontrivial background field. It is this background field 
which suppresses both the source and collision terms, to 
give 77/T 3 ~ |^| 2 at small I. 

We suggest that this is not an artifact of our approxi- 
mations. If deconfinement truly represents the ionization 
of color charge, then it is reasonable to expect that the 
propagation of all colored fields are suppressed at small 
£. In particular, while we have computed only to leading 
order in g 2 and l/log(l/g), one can show that within 
our ansatz, that 77/T 3 ~ £ 2 as I — > 0, order by order in 
perturbation theory at N c = 00 [29]. The suppression 
of r]/T 3 near T c is testable through simulations on the 
lattice [14]. 

If color is only partially ionized in the semi-QGP, then 
some color singlet states persist for a range of tempera- 
tures above T C! before they eventually disassociate into 
colored constituents. Similarly, with dynamical quarks 
colored fields contribute below T c , in the hadronic part 
of the semi-QGP. Including both effects is clearly chal- 
lenging; here we just speculate about the complete re- 
sult. Although hydrodynamics depends upon 77/s, with- 
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out a complete theory we can only discuss the behavior of 
i]/T 3 . Working up from low temperatures in the hadronic 
phase, r]/T 3 appears to decrease with increasing T, as in 
a liquid; e.g., for massless pions, to leading order in chiral 
perturbation theory [3, 9, 10]. Working down from high 
temperatures, from the complete to the semi-QGP, we 
have shown that the contribution of quarks and gluons 
to the shear viscosity decreases as rj/T 3 ~ \£\ 2 /a 2 (this 
assumes that higher order corrections in a s do not alter 
the leading behavior, at least qualitatively). Combining 
these two effects, we obtain a minimum for rj/T 3 near T c , 
as in non-relativistic systems [8, 15]. It seems implausible 
that the minimum is precisely at T c ; certainly in QCD, 
where there is no true phase transition [11, 13]. 

It is natural to suspect that heavy ion collisions at 
RHIC have probed some region in the semi-QGP. Since 
one needs a small value of the shear viscosity to fit the 
experimental data [1, 2, 3, 4, 5, 6, 7, 8, 9, 10], perhaps 
one is near T c . Heavy ion collisions at the LHC may 
probe temperatures which are significantly higher, pos- 
sibly well into the complete QGP. If so, then at small 
times collisions at the LHC create a system with large 
shear viscosity; as the system cools through T c , the shear 
viscosity then drops. Thus the semi-QGP predicts that 
at short times, the hydrodynamic behavior of heavy ion 
collisions at the LHC is qualitatively different from that 
at RHIC. This, and other phenomenological implications 
of the common suppression of colored fields at small £, 
will be developed separately [29]. 

In contrast, models of a strong QGP predict that hy- 
drodynamics at the LHC will be similar to that at RHIC, 
characterized by a small value of 77/s [3, 7, 10]. In par- 
ticular, while in JV = 4 supersymmetry the pressure/T 4 
is constant [16], several models have been proposed to fit 
the QCD pressure right down to T c [17]. Even so, in all 
of these models rj/ s is independent of temperature, and 
so remains small [17]. 

We eagerly await the experimental results for heavy 
ions from the LHC, which may be as unexpected and 
exciting as those from RHIC first were. 

This manuscript has been authorized under Contract 
No. DE-AC02-98CH10886 with the U. S. Department of 
Energy. 
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